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Abstract 


We analyze a system of linear algebraic equations whose solutions lead to a proof of 
a generalization of Boole’s formula. In particular, our approach provides an elementary 
and short alternative to Katsuura’s proof of this generalization. 
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1 Introduction 


Due to their numerous relations, binomial coefficients play an important role in various math- 
ematical fields, including enumerative combinatorics, statistics and number theory. In Boole’s 
classical book “Calculus of Finite Differences” [5], the following beautiful formula is given, 
which holds for 1 <m<neEN: 


Syey#(*)ee = a ae a) 


a ifm <n. 


It is known that the formula is related to Stirling’s partition numbers S(m, 7) (see, e.g.,[7]), 
which is given by the following equation 


S-y#(") k™ =n!-S(m,n) (m,n €N). 


k=0 


The relations between the formulas have implications for the study of the degrees of normed 
null-polynomials and the derivation of inequalities associated with the Smarandache function 
(see [10, 11]). 

The enduring interest in Boole’s formula has led to a variety of proof techniques being 
developed over the years. Gould [6] discussed its properties and called it Euler’s formula. 
In 2005, Anglani and Barile [2] introduced two proofs via methods from real analysis and 
combinatorics. Subsequently, Phoata [9] and Katsuura [8] provided new proofs and gave a 
generalization of Boole’s formula. In this note, we give a short and elementary proof of this 
formula which is based on a system of linear equations. More recently, Alzey and Chapman 
[1] have presented a novel proof, while Batir and Atpiar [3, 4] have independently developed 
two entirely new approaches to validating Boole’s formula. 
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2 Solutions of linear algebraic equations 


Let us consider the following case 


1 1 tee 1 Xo 0 
a a+b «::: a+nb Ly 0 
a2 (a+b)? =» (a+nb)?| |ro| — 0 |. (2) 
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where a,b can be any real numbers and the coefficient matrix V is a Vandermonde matrix. 
To solve this system, we first calculate the determinant of V. Since V is a Vandermonde 
matrix, its determinant can be computed as follows: 
: : n(n+1) 
det(V) = |] ((a+%b)- (a+ jb)) =nl-(n-1)!-- be. 
0<j<i<n 


We then proceed to define matrix V;, as follows: 


1 1 gas 1 0 1 re 1 
” a at+bo- at+t(k—1)b 0 at+(k+1)b --- a+tnb 
a” (a+b)” --- (a+(k—1)b)” b"-n! (a+(k4+1)b)” --- ae 


and we denote the matrix obtained by removing the (n+1)-th row and (&k+1)-th column from 
Vy, as Vey. The determinant of Vz; is computed by first applying Laplace’s expansion along 
the (k + 1)-th column, which yields the determinant of the submatrix V;;. This submatrix is 
also a Vandermonde matrix, and its determinant can be computed using the way previously 
showed. 


det (Ve) = (—1)"~* - b” - n! - det (Vir) 
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It is sufficient for us to use Cramer’s rule to obtain an explicit expression for xz 


= S0¥6) _ Cat = ay#(2), 


wk det(V)  -kl-(n—k)L k 


Reading the equations in (2) one by one we obtain the following result. 


Theorem. (Generalization of Boole’s formula) For any real numbers a and b, and for 1 < 
m<néEN, we have 


_ n m_ }(-1)?-6"-n! ifm=n, 
-() (a+ oa) xt 3) 


rh ifm<n. 


For the special case a = 0 and 6 = 1, the result in the Theorem implies Boole’s formula 
(1): 
Perhaps similar approaches from linear algebra can also be used to generalize other com- 
binatorial identities. 
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